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The Higgs boson production cross section at pp and pp colliders is calculated in QCD at next- 
to-next-to- leading order (NNLO). We find that the perturbative expansion of the production cross 
section is well behaved and that scale dependence is reduced relative to the NLO result. These 
findings give us confidence in the reliability of the prediction. We also report an error in the NNLO 
correction to Drell-Yan production. 
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INTRODUCTION 

Gluon fusion will be the most important production 
channel for Higgs discovery at the LHC. The Higgs boson 
should manifest itself in the reaction — > H{—> 77) +X, 
where a signal should emerge on top of a very smooth, 
measurable 77 background. At the Tevatron, the focus 
for Higgs discovery is in associated production modes like 
W/Z + H and tt + H. In a mass window around the WW 
threshold, however, gluon fusion is important. 

NNLO corrections to the process gg ^ H have been 
evaluated recently in the heavy top limit and the ap- 
proximation of soft gluon radiation ||, ^, where the 
partonic center-of-mass energy is close to the Higgs mass, 
Mfj/s = X ^ I. If we write the partonic cross section as 
an expansion in (1 — x), it has the following form: 
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where the [ ]+ terms are "+" -distributions defined in the 
usual way (see, e.g. Ref. Refs. [0, ||, ^ contain the 
coefficients a'-"-' and t^"'' up to n = 2 of this expansion. 
However, as anticipated in Ref. [Q, these contributions 
are not sufficient to arrive at a reliable prediction for 
the total cross section. Using resummation techniques, 
the authors of Ref. Q evaluated the coefficient at 
NNLO. It was included in the final results of Refs. |j. 
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However, the unknown sub-leading terms Cq/ with ? < 2, 

were treated in different ways by Refs. and leading 

to significant deviations in the numerical results. It is the 

purpose of the current letter to report on the analytical 
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evaluation of the coefficients Cj^, with k — 0, ... ,3 and 
Z > 0. In other words, we compute the partonic cross 
section for Higgs production in terms of an expansion 
around the soft limit. We find that the series converges 
very well and conclude that our final results are equiva- 
lent to a calculation of the cross section in closed analytic 
form. We therefore resolve the ambiguities of Refs. ||, ||] 
and provide a realistic prediction for the Higgs produc- 
tion cross section in pp and pp collisions. 

In checking our methods, we found an error in the 
NNLO Drell-Yan calculation of Ref. The correct re- 
sult is given at the end of the next section. 



THE CALCULATION 

In the following we will assume all quark masses to 
vanish, except for the top quark mass, and neglect all 
electro- weak couplings. In this limit, the Higgs boson 
can couple to gluons only via a top quark loop. This 
coupling can be approximated by an effective Lagrangian 
corresponding to the limit mt 00, which is valid for a 
large range of Mh, including the currently favored region 
between 100 and 200 GeV. The effective Lagrangian is 
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where G"^ is the gluon field strength tensor, H is the 
Higgs field, V « 246 GeV is the vacuum expectation value 
of the Higgs field and Ci{as) is the Wilson coefficient. 
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Renormalization of this Lagrangian has been discussed 
in Ref. for example, and shaU not be repeated here. 
In the MS scheme, the coefficient function Ci(as) reads, 
up to the order required here |||, |^ : 
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where l^t — ^'^■{fJ.'ji/M^). jXR is the renormalization scale 
and Mt is the on-shell top quark mass, a, = ai^\fij^) 
is the MS renormalized QCD coupling constant for five 
active flavors, and n/ is the number of massless flavors. 
In our numerical results, we always set = 5. 

The Feynman diagrams to be evaluated for hadronic 
collisions at NNLO are: (i) two-loop virtual diagrams for 
gg —f H; (ii) one-loop single real emission diagrams for 
gg Hg, gq Hq, and qq Hg; (iii) tree-level dou- 
ble real emission diagrams for gg — s- Hgg, gg — > Hqq, 
gq Hgq, qq -> Hqq, qq -> Hgg, and qq Hqq. The 

coefficients a*^^) and in Eq. (|l]) are determined by the 
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gg sub-process only, while the C;^, receive contributions 
from all sub-processes. 

For the single real emission diagrams (ii), the full an- 
alytical result for general values of x has been evaluated 
and will be published elsewhere. It can be expanded triv- 
ially in terms of (l — x). In order to obtain this expansion 
for the double real emission contribution (iii), we evalu- 
ated the squared amplitude and expressed the invariants 
of incoming and outgoing momenta, as well as the phase 
space measure in terms of two scattering angles and the 
dimensionless variables x, y, z, defined by M 
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Pi, P2, and ph are the momenta of the incoming partons 
and the (outgoing) Higgs boson, respectively. Modulo 
powers that vanish as d ^ 4 (d is the space-time dimen- 
sion), the resulting expression is expanded as a Laurent 
series in (1 — a;). The leading terms in this expansion are 
of order (1 — x)~^ and give rise to the purely soft contri- 
bution obtained in Refs. ||, ||]. Here we also keep higher 
orders in this expansion, (l — x)', / > 0. Aside from 
a few extra algebraic manipulations of hypergeometric 
functions, this expansion procedure allows us to perform 
the phase space integration along the lines of Ref. 0. 
Details of the calculations will be presented elsewhere. 

There are a number of checks that can be performed on 
our result. One is to see that all poles in d — 4 cancel. We 
have explicitly verified this cancellation through order 
(1 — xY^ . Since we have computed single real emission 
and the mass factorization counterterms in closed form 



(as opposed to an expansion in (l — x)), we can also 
obtain the pole terms for double real emission in closed 
form by demanding that the poles cancel. This allows 
us to obtain in closed form all finite terms in the cross 
section that are linked to the poles. These include all 
terms proportional to ln"(l — x) {n = 1,2,3) and all 
explicitly scale-dependent terms. 

As another check on our approach we applied it to the 
cross section for the Drell-Yan process at NNLO, where 
the full x-dependence is known in analytical form 
Detailed comparison with unpublished intermediate re- 
sults ijl^ shows that our expansion of the tree-level dou- 
ble real emission terms is in complete agreement with 
the corresponding expansion of the exact calculation. 
However, we find differences in the one-loop single real 
emission terms which we also have computed exactly. We 
conclude that the NNLO result for the Drell-Yan process 
in Ref. [p| is incorrect and that the correct result is 
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The numerical effect of these corrections is rather small 
and shall be investigated in more detail elsewhere. 

PARTONIC RESULTS 



We now present the result for the partonic Higgs pro- 
duction cross sections at NNLO. We define 
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The lower order terms, ct,-?^ and M]"^ , are given in 
Refs. M, pi. If we split the second order terms into "soft" 



and "hard" pieces, 

"ii — "w"3g "gg ^ "ii ' 
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the soft pieces are given in Eq. (25) of Ref. ||, while the 



hard pieces, ^l"^''^ (to order (1 — x)^) are: 
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For the sake of brevity, we have suppressed explicitly 
scale dependent terms by setting — ^J■R — Mh (they 
can be readily reconstructed using scale invariance) and 
displayed terms only to order (1 — x)^. Terms to order 
(1 — a:)^ dominate the corrections (see Fig. (^)), but we in- 
clude terms to order (1 — a;)^^ for all sub-processes in our 
numerical analysis. The labels "NS" and "S" in Eq. (|l^) 
denote the flavor non-singlet and singlet quark contribu- 
tions, respectively. The four contributions are equal only 
to order (1 — a;)^; their expansions differ at higher orders 
of (1 — x) (except that al^^g = (Tg^^g^ exactly). We note 
in passing that our explicit calculation confirms the value 
for the coefRcient cj^^ for the gluon-gluon subprocess de- 
rived in Ref. [|| . 

HADRONIC RESULTS 

The hadronic cross section a is related to the partonic 
cross section through a convolution with the parton dis- 



tribution functions. It has been argued that conver- 
gence is improved by pulling out a factor of x from dij 
before expanding in (1 — x). We indeed observe a more 
stable behavior at low orders of (1 — x) and will adopt 
this prescription in what follows. Beyond fifth order, 
however, it is irrelevant which is used. 

In Fig. (|l|), we show the cross section at LO, NLO and 
NNLO. At each order, we use the corresponding MRST 
parton distribution set ||l|] ||ll|, |l|l. The NNLO distri- 
butions are based upon approximations of the three-loop 
splitting functions |l3| . Studies using other parton distri- 
butions, including the NNLO distributions of Alekhin |l^ 
will be presented elsewhere. 
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FIG. 1: LO {dotted), NLO (dashed) and NNLO (solid:) cross 
sections for Higgs production at the LHC (hf = Hr = Mh)- 
In each case, we weight the cross section by the ratio of the 
LO cross section in the full theory (Mt — 175 GeV) to the LO 
cross section in the effective theory (Eq. (^). 

We next look at the quality of the expansion that we 
use for the evaluation of the NNLO corrections. Fig. (||) 
shows the NNLO if-factor (i^NNLO ^ ^nnlo/^lO) ^j. 
the LHC starting from the purely soft limit oc (1 — x)~^ 
and adding successively higher orders in the expansion in 



(1 — a:) up to order (1 
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Clearly, the convergence is 



very good: beyond order (1 — a:)^, the curves differ by less 
than 1%. Observe that the purely soft contribution un- 
derestimates the true result by about 10-15%, while the 
next term in the expansion, c>c (1 — a:)°, overestimates it 
by about 5%. Note that the approximation up to (1 — a:)" 
is not the same as the "soft -|-sl" -result of Ref. or 
the "SVC"-resuh of Ref. |], since these include only the 
ln''(l — x) terms at that order. 

We next consider the renormalization scale ((J-r) and 
factorization scale [hf] dependence of the if- factors. At 
the LHC, we observe that the fip and dependence has 
the opposite sign. In order to arrive at a conservative 
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FIG. 2: JsT-factor for Higgs production at the LHC. Each hne 
corresponds to a different order in the expansion in (1 — a;). 
The renormalization and factorization scales are set to Mh- 



estimate of the scale dependence, we display two curves 
corresponding to the values (/i^j, /iF) = {2Mh, Mh/2) 
and {Mh/2,2Mh) (see Fig. (|)). The scale dependence 
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FIG. 3: Scale dependence at the LHC. The lower curve of 
each pair corresponds to — 2Mh, IJ-f = Mh/2, the upper 
to /i_R = Mh/2, hf ~ 2Mh- The factor is computed with 
respect to the LO cross section at = fj,F = Mh- 

is reduced when going from NLO to NNLO and, in con- 
trast to the results in Ref. [Q, the perturbative series up 
to NNLO appears to be well behaved. The reason is that 
both the newly calculated contributions from hard ra- 
diation and the effect of the previously unavailable set 
of NNLO parton distribution functions reduce the NNLO 
cross section. Detailed studies of the individual effects 



will be presented in a forthcoming paper. 

Fig. ^ shows the results for the Tevatron at a center-of- 
mass energy of ^/s = 2 TeV. Here the dependence on 
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FIG. 4: Scale dependence for Tevatron Run II. The lower 
curve of each pair corresponds to fiR = fiF = ^Mh, the upper 
to fj,R = ^iF = Mh/2. 

and has the same sign, so we set /i^ = fip = ^ and 
vary fi between Mh/2 and 2Mh- The if-factor is larger 
than for the LHC, but the perturbative convergence and 
the scale dependence are satisfactory. 



CONCLUSIONS 

We have computed the NNLO corrections to inclusive 
Higgs production at hadron colliders. We find reasonable 
perturbative convergence and reduced scale dependence. 
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